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Abstract 

We have derived effective gravitational field equations on a lower dimensional hypersur¬ 
face (known as a brane), placed in a higher dimensional bulk spacetime for both Einstein 
and /(TT.) gravity theories. We have started our analysis on n-dimensional bulk from which 
the effective field equations on a (n — l)-dimensional brane has been obtained by imposing 
Zi symmetry. Subsequently, we have arrived at the effective equations in (n — 2)-dimensions 
starting from the effective equations for (n — 1) dimensional brane. This analysis has been 
carried forward and is used to obtain the effective field equations in (n — m)-dimensional 
brane, embedded in a n-dimensional bulk. Having obtained the effective field equations in 
Einstein gravity, we have subsequently generalized the effective field equation in (n — m)- 
dimensional brane which is embedded the n-dimensional bulk spacetime endowed with f{TZ) 
gravity. We have also presented applications of our results in the context of Einstein and 
f{TZ) gravity. In both the cases we have discussed vacuum static spherically symmetric 
solutions as well as solutions in cosmological context. Implications are also discussed. 


1 Introduction 

One of the fundamental problems in theoretical physics amounts to unifying all known inter¬ 
actions. A consistent unification of gravity with the other fundamental interactions has been 
an active area of research over many years. One of the leading contenders along this direction 
is superstring theory. In superstring theory the universe is assumed to be 11-dimensional of 
which seven are compactified leaving four non-compact dimensions, which we observe. Among 
various candidates, the 10-dimensional Eg, x Eg heterotic string theory is a strong candidate as 
the theory is able to accommodate the standard model within it. Recently it has been shown 
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that the 10-dimensional Eg, x Eg heterotic string theory is related to the 11-dimensional M 
theory (with an orbifold symmetry S^lZ 2 ). This also confines the standard model particles to 
standard non-compact 4-dimensional spacetime while gravity can probe the higher dimensional 
spacetime [1-3]. 

A simplistic model which captures most of the key notions is a 5-dimensional model in which 
matter fields are confined to 4-dimensional spacetime while gravity exists in all the five spacetime 
dimensions. The first step along this direction with this simplified picture presented above 
was taken by Randall and Sundrum by introducing two 4-dimensional sub-manifolds (known 
as brane) in a 5-dimensional anti-de Sitter bulk [4]. This arrangement of two positive and 
negative tension branes separated by a finite distance (known as radion field [5-7]) can address 
the hierarchy problem by lowering the Higgs mass to electroweak scale on the visible brane. 
Subsequently, they have introduced another model, which consists of a positive tension brane 
with an infinite extra dimension [ 8 ]. Due to the presence of an extra spacetime dimension it is 
expected that there would be deviation of various physical results from Einstein gravity, which 
would become more significant in the high energy limit. The brane world model subsequently, 
was applied in the context of particle phenomenology [9-12], black hole physics [13-17] and 
cosmology [7,18-21]. 

A natural extension of the Randall-Sundrum (RS) scenario to models with more than one 
warped extra dimension have been proposed, where several independent jZ^ orbifold is con¬ 
sidered along with the 4-dimensional non-compact manifold M 4 [22-26] . These multiply warped 
models are mainly motivated from various considerations, which include: radion stabilization, 
fermion mass hierarchy in Standard Model etc. There have been significant number of follow-up 
works in these multiply warped models which involves: inflation, cosmic acceleration, matter 
field localization, Kaluza-Klein modes of graviton, gauge and scalar fields [27-32]. However all 
these results depend on certain assumptions of bulk metric elements and bulk Einstein equation 
under some simplified conditions, for example in presence of homogeneity and isotropy. 

There exist an elegant way in which effective field equation for gravity on a lower dimensional 
spacetime can be obtained. This involves use of geometric quantities like metric, curvature 
induced on the brane from bulk. Using these induced geometric quantities it is possible to 
arrive at a relation between the bulk curvature tensor and the brane curvature tensor, known as 
Gauss-Codazzi equation [33,34]. From various contractions of the Gauss-Codazzi equation we 
can determine the effective field equation on the brane by relating bulk Einstein tensor to that 
on the brane. 

This approach was applied in order to obtain the gravitational field equation on the brane 
from a 5-dimensional bulk in [35] . The effective equation contains non-local terms inherited from 
the bulk. Subsequently starting from the effective field equation derived in [35] a class of vacuum 
solutions has been obtained in [16], which has been further generalized in [36]. 

It is strongly believed that general relativity is only a low energy approximation of some 
underlying fundamental high energy theory [37, 38] . This suggests that Einstein-Hilbert action 
in high energy limit, should be modified by introduction of higher curvature terms. Such an 
alternative theory is the f{TZ) theory, where TZ is the bulk Ricci scalar. This theory has been 
studied extensively in the context of solar system tests, inflationary paradigm, late time cosmic 
acceleration along with possibility of detecting gravitational waves [39-43] . Moreover two-brane 
models, due to the presence of JiTZ) term in the bulk has also been a topic of discussion recently 
in collider physics [44-46] relating to the absence of graviton Kaluza-Klein Modes in LHG. 

The above approach of obtaining effective gravitational field equation is usually confined to 
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Einstein gravity on the 3-brane alone. This subsequently has been generalized to obtain effective 
field equation on the 3-brane when the bulk contains f(TZ) gravity. This effective field equations 
has been used in [47] to obtain vacuum solution in f{TZ) bulk. 

In this work, we generalize the above setup by working in an n-dimensional bulk and then 
obtaining effective equation on an (n — m)-dimensional brane where m can take any values. 
This seems natural from a string theory view point, where the universe is 11-dimensional, while 
our universe is a 3-brane with 4-spacetime dimensions. Following this argument we have derived 
effective gravitational field equation on a (n—m)-dimensional brane embedded in a n-dimensional 
bulk. We start with a n-dimensional bulk spacetime and then derive the effective field equation 
on a (n — l)-dimensional brane, which can then be used as the starting point to derive effective 
field equation on a (n — 2)-dimensional brane. Having obtained the effective field equation on 
the (n — 2)-dimensional brane one can easily recognize the pattern of obtaining the effective 
field equation in any arbitrary lower dimensional surfaces. Thus following the road from n- 
dimensional bulk to (n — 2)-dimensional brane, we can obtain the effective field equation on 
(n — m) dimensional brane as well. The above analysis can be easily generalized to f{TZ) gravity 
model in the bulk following the analysis presented in [47]. 

The paper is organized as follows: In Sec. 2 we have elaborated the derivation of effective field 
equation for Einstein gravity. First, in Sec. 2.1 we present the derivation of effective field equation 
in (n — l)-dimensional brane, which is subsequently generalized in Sec. 2.2 to (n — 2)-dimension. 
Following the pattern we obtain the effective gravitational field equation in (n — m)-dimensional 
brane in Sec. 2.3. This result has been generalized for f{TV) gravity model in Sec. 4. The effective 
equations so derived have been applied to obtain static spherically symmetric vacuum solutions 
and cosmological solutions along with their compatibility with standard cosmological evolutions. 
The results in the context of Einstein gravity have been discussed in Sec. 3 while those for f{TZ) 
gravity is discussed in Sec. 5. Finally we conclude with a discussion on the results obtained. 

Throughout the paper, we maintain the following convention: all Latin letters a,b,... run 
over all the n spacetime indices, Greek letters a,/r,. .. run over the (n — 1) spacetime indices. 
Finally capitalized Latin letters H, H,... run over the (n — 2) spacetime indices. 

2 Effective Field Equation In Einstein Gravity: Background 
and Formulation 

Hamiltonian formulation in general relativity comes up with the (1-1-3) splitting of the spacetime, 
where a time direction is singled out. For this purpose we introduce a time coordinate t = t{x^) 
and foliate the spacetime with such t = constant surfaces. Introducing a vector field t = d/dt 
we can move out of the surface by moving along the integral curves for t. This involves a change 
in the t coordinate introducing the shift function N and change in the orientation captured 
by the lapse functions 7V“. Then TV, yield 4 components of the Einstein tensor, while the 
remaining six components are being determined by induced metric on the surface hap. The 
induced metric can be obtained from the bulk metric gab in two ways. First, we can introduce 
normalized normals to the surface given by: Ua = —N'Vat and then define the induced metric as 
j^ab _ gab_^a^b formulation does not work for null surfaces, for which a complete treatment 
has been provided in [48]). Secondly, we can introduce coordinates on the t = constant 
hypersurface on which we can introduce surface tetrads e“ = dx°‘/dy°^ and obtain induced 
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metric as: hap = gabSas’p- With the induced metric we can define covariant derivative on the 
surface and then using vectors on the t = constant surface we can introduce curvature tensor 
on the spatial hypersurface. The relation between 4-dimensional and 3-dimensional curvature 
tensor is obtained through Gauss-Codazzi equation. This can be manipulated in various ways 
in order to arrive at the effective Einstein’s equation on the t = constant hypersurface. This 
procedure can be extended easily for timelike surfaces (i.e. spacelike normals) as well. 


2.1 Effective Field Equation on (n-l)-diniensional Brane 


Having described the general setup let us now consider the effective field equation on a brane 
embedded in a higher dimensional bulk, which is the main motivation of this work. The procedure 
sketched above can then be applied for a generic situation where a single (n — l)-dimensional 
brane is embedded in a n-dimensional bulk spacetime with Einstein gravity. Eor this set up the 
effective gravitational field equation on the (n — 1) dimensional brane can be written as (see 
App. A.l): 




I n — 3 
‘n- 2 


Tab^a^^p + ^aP fTabn°‘n'^ — 


^ab 


n — I 


-T 


-eE, 




-b e 


KK^p - - -h^p {K^ - 


( 1 ) 


The above expression connects bulk energy momentum tensor and the extrinsic curvature to the 
(n— 1) dimensional Einstein tensor. Given Tab, along with surface contribution we can solve the 
above equation to derive (n — 1) dimensional solution. 

The above effective equation is applicable to both spacelike and timelike surfaces. However 
for our purpose, we will need only timelike surfaces, which are the branes. Thus the basic 
assumption to be followed in this section is that the n dimensional spacetime can be projected 
on a (n — 1) dimensional brane with spacelike normal. The line element then takes the form: 

ds^ = dx^ + hapdx°'dx^ ( 2 ) 


The brane is located at y = 0 on this n dimensional spacetime. Using the standard assumption 
that normal matter can exist only on the brane, the bulk energy momentum tensor can be 
written as: 

Tab = -Kgab + SapCae^S (y) (3) 

where, Sap represents the brane energy-momentum tensor and A„ is the bulk cosmological 
constant. As the surface y = 0 contain surface energy momentum tensor, the extrinsic curvature 
is not continuous. The discontinuity is related to the surface stress energy tensor through the 
following relation: 

Sap = — {[Kap] — haP [K]) (4) 

where and [K\ denote jump in Kap and K across the y = 0 surface. Thus we obtain the 

jump in the extrinsic curvature which after imposing Z 2 symmetry, has the expression 

^aP ~ ~^aP ~ 2 '^n (^SaP — ^ _ ^hapS^ (5) 
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From the extrinsic curvature we can contract both the indices to obtain a scalar. This scalar 
K corresponding to each side of the x = 0 hypersurface is represented by K^. Therefore from 
Eq. (5) it has the expression; = —K~ = — (^K^S/2{n — 2)). Using both the expressions for 
Kajs and K we finally arrive at the following identities: 
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(9) 


In order to proceed further we need to decompose the brane stress-energy tensor in two parts. 
One of them involves the brane cosmological constant while the other provides matter stress- 
energy tensor on the brane. Implementation of this separation leads to, 


SaP — ^hap Tap 


( 10 ) 


where, a represents the brane tension. Next we substitute this decomposition of Sap into the 
above expressions involving extrinsic curvature. This ultimately leads to. 
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( 11 ) 


Eqs. (3) and (11) now need to be substituted in Eq. (1) in order to obtain the final form of the 
effective equation on the brane. After the substitution the effective Einstein equation on the 
brane takes the form 
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where we have introduced the second rank tensor Ila^ as: 

—^Vp + + 4(„_2)' 
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hoi0T 


( 12 ) 


(13) 


After rewriting Eq. (12) can also be presented in a more compact form, which can be expressed 
as: 
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where the energy momentum tensor has the expression: 


^Q/3 -^71 — 1^0/3 '^OL^ H” 


^a/3 


-En 


(15) 


In the above expression A„_i represents the (n — l)-dimensional cosmological constant which 
has the expression, 


An—1 


(n-3) kI 
(n- 1) kI_^ 



71—1 

8(n-2) 


2 2 


(16) 


This generalizes the expression for the effective cosmological constant on a lower dimensional 
hypersurface embedded in a higher dimensional bulk. For n = 5 we retrieve fine balancing 
relation of the Randall Sundrum single brane model [35] . 

Thus this formalism, with a negative n-dimensional cosmological constant can provide a 
solution to the cosmological constant problem by tuning the brane tension such that An-i = 0. 
In which case the brane tension can be obtained as: 

2 8(n — 2)|A„| 

Note that once the (n— l)-dimensional cosmological constant has been set to zero it will remain 
zero in all the lower dimensional branes. Also G„_i is the (n — 1) dimensional gravitational 
constant having the expression, 

kLi = 87rG„_i = (18) 

Having discussed the effective gravitational field equation in (n — l)-dimension, we will now move 
one more step by deriving effective field equation in (n — 2)-dimensional spacetime. 


2.2 Effective Field Equation on (n-2)-dimensional Brane 

Similar analysis can be performed while obtaining effective field equation on the (n — 2) dimen¬ 
sional brane starting from the (n — l)-dimensional spacetime. The field equation in the (n — 1) 
dimensional spacetime can be used which ultimately leads to the following expression for effective 
field equation in (u — 2) dimensions as (see App. A. 2): 
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AB 
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This expression represents the effective Einstein’s equation in an (n — 2) dimensional space- 
time. Our next task is to rewrite the above equation in terms of various components of energy 
momentum tensor, especially the extrinsic curvature components. 
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Let us now try to match our solutions to that of (n —2) dimensional brane energy momentum 
tensor. As usual we assume the following form of metric ansatz, 

= dx^ + hapdx°‘dx^ 

= dy^ + ha^e’^e^dQ^ + qAsdx^dx^ (20) 

where qAB is the induced metric on the (n — 2)-dimensional brane. Following the previous 
derivation of effective field equation here also we divide the energy momentum tensor, Tap in 
two parts, tap pertaining to (n — 1) dimensional brane and ^ab as energy momentum tensor in 
the (u — 2) dimensional brane. Thus we have the following relation, 

Fa/3 — tap ~\~ S ^AB^a ^/3 

= tap + S {() eaep {-T,qAB + I'Ab) (21) 

where E is the brane tension and i>ab is the brane energy momentum tensor. Then the tensor 
^ab can be obtained from the discontinuity of the extrinsic curvature from the relation, 

^AB = —— {[K,ab] — qAB [Aii]) (22) 

The jump in the extrinsic curvature by imposing Z 2 symmetry can be obtained in terms of ^ab 
as, 

^ab ~ ~^AB ~ 2 ^"“^ ^ _ ^Qab^ (23) 

Having obtained the extrinsic curvature, the scalar K, can be obtained by contraction with qAB- 
This has the expression: /C+ = —/C“ = — («:^_;^^/2(n — 3)). We therefore arrive at the following 
expressions for various combinations of ^ab as they appear in effective field equation: 

i^AB = {n — 2)Ti^qAB — vYiqAB — {n — 2)'Zi >ab + vvab (24) 

^ac^b = ^^Qab — ‘2-'ZiyAB + bac^b (25) 

^2 = 2)2^2 - 2(n-2 )Ez/ + z/2 (26) 

= {n- 2)Y? - 2 Sb + lyABi^^^ (27) 


This enables us to write an expression for the extrinsic curvature contractions in terms of ^ab, 
which can further be simplified in terms of i^ab following the above identities as: 
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It now just needs some rearrangement to write down the effective field equation on the (n — 2) 
dimensional brane. After inserting all relevant expressions and then rearranging the effective 
Einstein’s equation on this (n — 2) dimensional brane we obtain, 
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where, we have introduced, the following quantities, 

1 1 


IfcK/S - 

Tab = 


+ -^hapt^yt^'^ + 


4(n - 2) 


ttrvfl 


8(n-2) 


h^pt 


vac^b 


1 


:qABVCDV 


.CD 


1 


4(n — 3) 


VVAB 


8(n- 3) 


qABV 


(29) 

(30) 

(31) 


The effective equation can be written in conventional form after some algebraic manipulations 
such that: 


^Gab = k^-2Eab 

where energy momentum tensor Tab has the expression: 
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Here, we have introduced the following two objects for notational simplicity. 
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where A „_2 and k „_2 can be interpreted as the (n — 2)-diniensional cosmological constant and 
(n — 2) dimensional gravitational constant respectively. However for the case, in which tap = 0, 
i.e., the bulk energy momentum tensor vanishes, the above equation gets simplified to, 

= kI_^Tab (36) 


where Tab has the simplified form: 


Tab = —A„-29ab + vab + 


’'n-l 


Tab — 


-£ab 


n — 4 
n — 3 


1 


Eape‘Xe^^ + QAB i Eaps°-s^ 



(37) 


Thus the structure of effective held equation in in — 2)-dimensional brane for vanishing bulk 
energy momentum tensor is mostly identical to the (n — l)-dimensional counterpart except for 
the Weyl curvature parts. The structure of the effective equation should be quiet evident at this 
stage, which we are going to exploit later to obtain effective held equation in (n —m)-dimensional 
brane. 


2.3 Generalization to (n-m)-dimensional Brane 

Having discussed the effective gravitational held equation in (n — 2)-dimension starting from 
n-dimension through a two step decomposition, it is now time to consider the generalization of 
the effective held equation for gravity to (n — TO)-dimension. This can be obtained following the 
analogy of previous results. The hnal expression for effective held equation in (n — to) dimension 
is: 
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where the energy momentum tensor has the following expression: 


^AB 


i=2 


(n—m) 

-rriQAB 

(n-m) 

+ ^AB + 

^n—m+1 

^n—m 

^ n — i — 2 



An—i+l) 


) 



- AB 


71 ^AB 


(n—m) 




KZ 








n — 2 

n — 2 


(38) 




(39) 


9 




















with the following identifications: 
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However it may be observed that in the situation where there is only matter on the (n — m)- 
dimensional brane but not on the other higher dimensional branes, all the additional terms 
vanishes except the Weyl tensor part. If the bulk spacetime is assumed to be Anti de-Sitter 
(AdS) then the Weyl tensor can taken to be zero, in which case the effective field equation takes 
a simplified form: 


(n—n 


^Gab 


— A I _("-"») 

^''■n-mqAB ^AB 
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“72 ^ AB 


(44) 


Thus in this simple situation the effective field equation has contribution from (n—m) dimensional 
cosmological constant, energy momentum tensor on the brane and higher order correction terms 
originating from Tab- Having derived the effective field equation in Einstein gravity let us now 
consider the identical situation in f{TZ) gravity as well. 


3 Applications in General Relativity 

Having derived the effective field equation on an arbitrary lower dimensional hypersurface start¬ 
ing from a n-dimensional bulk, in this section we will consider two applications of our result. 
One of these applications would be for static, spherically symmetric configuration and the other 
in cosmology. As an illustration we will take the bulk to be six dimensional while the brane is 
obtained by two step decomposition and is four dimensional. 


3.1 Static Spherically Symmetric Vacuum Brane 

We start by considering the case for static, spherically symmetric spacetime with no matter 
source present. Thus brane energy momentum tensor vanishes however the electric part of the 
Weyl tensor survives. Hence the effective gravitational field equation on the four dimensional 
brane starting from the six dimensional bulk can be obtained from Eqs. (38) and (39) as. 


^‘^^Gab = —£ab 


2 

3 





(45) 


In arriving at the above equation the cosmological constant has been neglected, since we are 
interested in a local spacetime region. The electric part of the Weyl tensor in a static, spheri¬ 
cally symmetric spacetime can be written using two unknown functions depending on the radial 
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coordinate as, 
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where qab is the induced metric on the t = constant surface, such that u^qab = 0. The five 
dimensional electric part of Weyl tensor would have similar decomposition depending on two 
separate radial functions U{r) and P{r). On projecting this five dimensional tensor Eap on the 
four dimensional brane, the effective field equations reduce to 
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Now Kn and Kn-i are related through Eq. (18), such that, = (16/3 )k|(1/(t), k| = ( 6 /S)k 4 
and K 4 = SttG. Using these relations the above equation can be written in terms of the brane 
tensions, a for the five dimensional brane, E for the four dimensional brane and the Gravitational 
constant G as, 
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Now assuming a static, spherically symmetric metric ansatz of the form, 

ds"^ = + r^dn^ 

we obtain various components of the gravitational field equations to be 
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Hence the effect of multiple extra dimensions is manifested in the lower dimensional effective 
gravitational field equations through the appearance of extra terms U(r) and P{r). Field equa¬ 
tions have to be supplemented by the conservation equations, which will provide a differential 
relation connecting t/(r), P(r), U{r) and P{r). The conservation equations turn out to be, 
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Starting from Eq. (50) it is straightforward to obtain a solution for e~^. Since all the objects on 
the right hand side of Eq. (50) depends only on the radial coordinate we can integrate them to 
obtain the following solution for e~^ as, 
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47rGE 


dr r^U{r) + 
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Having obtained e~^ it is now pertinent to consider the other metric component e'^ as well. This 
involves all the extra contributions from dark radiation and dark pressure terms. In order to 
obtain solutions for these higher dimensional contributions, we need to impose certain relations 
among these higher dimensional objects. The most useful among them being the relations 
2U + P = 0 and 3U + P = 0. Then from Eq. (53) and Eq. (54) we readily obtain. 


P{r) 

Qir) 


-2U{r) = % P{r) = -3U{r) = 


Qo + 


3aPo 

2r 


PqP , 
9r3 ’ 
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47rGE’ 


h 

J.Q 
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32 TE 
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(56) 

(57) 


Eor these relations between dark radiation and dark pressure the static spherically symmetric 
solution corresponds to. 


2GM + Qo _ 3aPo _ ^ 
r 2r‘^ 9r‘^ 


(58) 


Note that this is an exact solution to the effective field equations, no assumptions or approxima¬ 
tions have been made. Hence the effect of extra spatial dimension appears through the charge 
term, i.e., r~^ term and the term. Thus near r = 0, 1/r^ term would dominate, which 
originates from the existence of more than one extra spatial dimensions. 

Having derived the static spherically symmetric solution there are two important aspects one 
can study — (a) stability of the solutions and (b) thermodynamics. We will briefly comment 
on these features. The stability of solution can be obtained by considering perturbations, which 
again are of three types, scalar, vector and tensor modes. These modes satisfy a set of decoupled 
wave equations with the structure. 


(□-7^r)* = o 


(59) 


Here □ represents d’Alembertian operator, /(r) = e" = e~^, $ stands for the perturbation 
modes and V is an appropriate potential associated with perturbation $. Eor stability, the 
potentials should be positive. If we assume that the corrections to the Schwarzschild solution is 
small, i.e., a and /3 are small compared to M, then the above solution is stable for all the modes. 

Let us now discuss the case of black hole thermodynamics. For the horizon location we have 
the following equation. 


2GM Qq 
rh 


3aPo PoP 


(60) 
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The mass of the black hole (or, equivalently the internal energy) can be obtained readily from 
the above equation as, 


M{rh) = y 


Qo 3 q;Po PqP 
2 Art 6rl 


(61) 


Now entropy can be obtained from holography, i.e., S = AjA, such that, S = r\. Thus internal 
energy can be obtained in terms of entropy density from Eq. (61) such that. 


M{S) 


Vs Qo 3aPo PqP 

2 2 4v^ 65'3/2 


Then the black hole temperature can be obtained in terms of entropy as, 



1 3qPo PoP 

aVS 8 ^ 3/2 455/2 


which immediately leads to the following expression for specific heat as. 


Cv=T 



1 I 3aPn I Po0 

AVS 4S5/2 

1 I 9aPo I 5Po/3 

^5372 -I- 16S5/2 -I- ^5772 


(62) 


(63) 


(64) 


Thus the specific heat of the black hole is negative, as it should and there is no phase transition 
possible, since a and j3 are both positive quantities. Later we will show that inclusion of f(TZ) 
gravity can alter the situation drastically, which we will discuss in the next section. 

Having derived the static spherically symmetric solution and its various properties, let us 
now discuss some novel features of this solution 


• Imposing the following restrictions 2U(r) + P{r) = 0 and P(r) + 3C7(r) =0, we have shown 
that the metric elements satisfy e~^ = e'^. Thus the line element derived in Eq. (58) 
depends explicitly on the form of U{r) and U(r). Thus the geometry of the spacetime is 
uniquely determined only when knowledge about multiple dimensions are obtained through 
U{r). Also it turns out that at small r the effect from U{r) dominates. 

• The solution is also stable under perturbations. The potentials associated with respective 
wave equations are positive for scalar, vector and tensor modes. Thus the black hole 
solution, will not change its character under perturbations. 

• Finally, the specific heat associated with the black hole solution is negative. This shows 
that the black hole is thermodynamically unstable. Since all the additional parameters are 
positive there is no phase transition in this black hole system. 

Hence the fact that spacetime has more than one extra dimensions is being reflected on the 
right hand side of the effective gravitational field equations. These extra terms contribute to the 
energy momentum tensor and modifies the metric components from just being 1 — (2GM/r). 
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3.2 Cosmological Spacetime 

After discussing the static spherically symmetric configuration it is natural to turn around and 
consider gravity theory on a lower dimensional hypersurface in a cosmological context. As 
already mentioned we will start with 6-dimensional bulk spacetime and then obtain cosmological 
solutions from the effective field equation on a 4-dimensional brane. In order to obtain analytical 
solution, we assume Eat = 0 i.e. contribution from electric part of Weyl tensor vanishes. Then 
effective field equation in 4-dimension takes the following form: 

Gab = ~ Ba^a<1ab + i^aTab + 

+ + QAB |c/3s“s'^ - 

-f {^apeAes + Qab |na;3s“s^ - (65) 

In this expression along with brane energy momentum tensor tab and its higher order term 
two additional contribution from 5-dimensional matter energy-momentum tensor is present. 
We assume that the brane is filled with perfect fluid with the following form for the energy- 
momentum tensor on the 4-dimension as: 

Tg = diag {-p,p, p, p) (66) 

where p represents energy density of the perfect fluid and p yields its pressure. The energy- 
momentum tensor in 5-dimension is taken to be pressure free dust, such that its energy momen¬ 
tum tensor has the simple form: 


= diag (-P 5 ,0,0,0,0) (67) 

where the subscript ‘5’ indicates that it is from 5-dimensional brane. Using which we arrive at 
the following expression for the higher order components 

Hit = n = 

Having obtained all the components of the matter energy momentum tensor, the time-time and 
spatial components of effective equations lead to (we will assume that the spacetime is spatially 
flat, i.e., k = 0): 


377 ^ = kIK^ + kIp + + ^4p5 ( 69 ) 

+ = - kIp - {P^ + 2PP) - (70) 

Having obtained the Friedmann equations, let us explore the various phases of the universe. To 
get a handle, let us write the above expression in terms of relative abundances fli = Pifi/pc- 
Here pi^ stands for the present energy density of the ith component and pc is the current total 
energy density of the universe. The energy density can be expanded in terms of baryonic and 
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relativistic matter fields. We will assume that matter in higher dimension to be non-relativistic. 
Thus the Hubble parameter turns out to be, 


= Hi 


Ha + (Hs 


05 ) 0 ’+ n«p 


\ a 


(71) 


where, H 5 = (K|/ 167 rG£)(p 5 _o/pc), with I being a characteristic length associated with the extra 
dimension. Note that the above Hubble parameter can reproduce the standard KCDM model. 
For at late times Ha would dominate leading to accelerated expansion of the universe. At early 
times, the 1/a^ term would dominate, leading to radiation dominated universe. The intermediate 
region is dominated by non-relativistic matter such that ^ a~^. However the surprising fact 
about this model is that in the matter sector along with the baryonic matter, due to extra 
dimension we have additional contribution, which can act as dark matter. In standard KCDM 
cosmology, the dark matter has to be included in the model by hand. However in this case it 
appears in the Friedmann universe naturally, we need not have to invoke it by hand. Using 
the epoch of matter, radiation equality it is possible to get a numerical estimate about five¬ 
dimensional gravitational constant as, 

! + = 3.9x10* (0.022 + ^'^) (72) 

where h = 0.72, I stands for the length of the extra dimension. Thus knowing the redshift 
Zeq, we can determine the value of the higher dimensional gravitational constant Thus 
spacetime having more than one extra spatial dimensions can naturally provide additional matter 
components in the right hand side of the Friedmann equations. To have inflation we need to 
include an additional scalar field on top of this alike the standard KCDM model. Thus the 
effective field equations are capable to explain the phases of the universe. 

Writing k^K^ = Aes and for the bulk matter we have ps = (p 5 _o/^)(l/a^), where ^ rep¬ 
resents the finite length of the 5-dimensional brane coordinate. Assuming that the matter in 
the 4-dimensional brane is also pressure-less, we have the following expression for the Hubble 
parameter: 




Aeff , SttC Po ^ 1 4 Pq ^ 1 2 P 5 .O 1 , 1 4P5,0 1 

3 3 a3 e"*® £ a3 64"'® £2 ^6 


(73) 


where po represents the value of energy density in the present epoch and p 5 ,o represents the 
same for 5-dimensional matter. The above expression for Hubble parameter can be simplified 
and arranged properly leading to: 
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Aeff 

“ 3 “ 


SttG 


-Po 


1 2P5,0 


1 


36 


4 2 
^ 5/^0 


64 



1 


(74) 


Then following the techniques developed in [49] we can introduce three constants, namely: 


Ca = —, 


SttG 

— o 7*0 


1 2^. ^ _J _^4 2 , J_ 4 ^ 27.3 

6 ^^ £ ’ 36 ^^°+ 64 "'® £2 ( 75 ) 
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such that we obtain the following solution to the scale factor as: 


exp 


3\/^ (t — to) 


2VCly'CAa^ + Cpa^ + Csq + 2CAa^ + Cp 


2'J Ca \J Ca + Cp + Csq + 2Ca + Cp 


Ca > 0 (76) 


3\/-Ca (to -t) = sin 
Ca<0 


-1 


2 Cao^ + Cp 

I Cl - 4CaC's„ 


C^p > 4CAC'sq 


— sin 


-1 


2Ca + Cp 


I Cl - 4C'ACsq 


(77) 


Though the solutions have a complicated appearance, we can impose certain conditions under 
which the solutions simplify substantially, yielding clearer physical insight. For that purpose we 
consider the situation Aes = 0. This would be true if we assume the 5-dimensional cosmological 
constant to be negative, such that: A 5 = —(k|E^)/ 6 . Under this condition the Hubble parameter 
from Eq. (74) turns out to be. 




sttg r i4 P5 
3 ^ IGttG p 


2 E ^ 16k|p2 ) 


(78) 


If we assume that p represents energy density of non-relativistic matter we arrive at the following 
differential equation 


2 ^ [ SttGpo ^P^\ J_ 47rGp§ / J_ 

\ 3 6 £ J 3E y 16/^5 

leading to the following solution for the scale factor: 


a 


3 


| 6 ,G« + 3 A^is|t^ + 


127rGpg / 9 pi,o \ 


t 


(79) 


(80) 


It is clear from the above expression that the universe undergoes a transition in the expansion 
rate. In the standard scenario there are three phases, the radiation dominated, the matter 
dominated and the accelerated expansion. If we treat the additional correction terms due to 
existence of extra spatial dimensions with importance then even within matter dominated epoch 
we will have transition. The time scale when this happened depends on whether the brane 
matter energy dominates over the bulk matter or not. Below we present the result for timescale 
of transition for both the situations: 


^ = —— = --A"^- 

SttGE EkI 3 ® ’ 

t (1 - 

V kIps.o ) 

Thus at early universe (within the matter dominated epoch) we have a high energy regime, where 
a ~ while at late time low energy regime the scale factor variation with time modifies to 

a ^ which is the standard evolution of the matter field. Note that the same should happen 

in radiation dominated epoch as well. Hence the scale factor at later times will behave as it 
should in a standard ACDM model. In a nutshell, the effect of extra dimensions on cosmological 
scenario can be seen from the following results: 


{P5,o/^ “C Po) 

(81) 

{P5,o/^ > Po) 

(82) 
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• It is clear from the expression for Hubble parameter that the 5-dimensional matter adds to 
the 4-dimensional one and enhances the total non-relativistic matter content of the universe 
(4-dimensional brane we live in). Hence matter fields present in the extra dimension can 
act as a viable source of dark matter. The most important part of the above analysis 
corresponds to the fact that the dark matter contribution appears naturally, we need not 
have to add them by hand (a related issue has also been addressed in [50]). 

• The existence of more than one extra spatial dimension introduces additional correction 
terms in the Friedmann equations. This results in transition within both matter dominated 
and radiation dominated epoch. This in principle can affect the structure formation in the 
universe and hence can be used to constrain various higher dimensional parameters. 


Thus with matter present in the 5-dimensional brane we can have standard cosmology with 
proper scaling of the scale factor with time. However the matter content and transition time 
scale will depend on the content in the extra dimensions. 

The cosmological solutions in brane models have interesting holographic interpretations. It 
was shown in [54] that the usual four dimensional cosmological equations can be re-written in 
an equivalent form of the entropy of two dimensional conformal field theory (CFT). In our case 
the Friedman equations are given by Eqs. (69) and (70) respectively. Now defining the Hawking, 
Bekenstein and Bekenstein-Hawking entropies as, 


{n-l)HV 
~ 4G ’ 


Sbh 


{n-l)V 

4Ga 


Ss = —E 

n 


(83) 


where E = pV is the total energy within a volume V and pm oc a“". Also the total entropy 
should be bounded by the Hubble entropy Sh, which is the entropy of the black hole with the 
radius of the Hubble size. Using these definitions we can rewrite the Friedman equations in the 
form, 


Sh = -V Ebh (2A — KEbh)', KEbh = n{E + pV — ThSh) (84) 

n 

with energy of the black hole being Ebh = n{n — l)U/87rGa^ and Hawking temperature being 
Th = —H/2ttE[. This is exactly the form one would have from the CV entropy relation. It 
explicitly shows the connection between four dimensional classical gravitational physics and 2- 
dimensional quantum 2-dimensional CFT [55] . This two dimensional CFT entropy representation 
of cosmological equations is called generalized Cardy or Cardy-Verlinde (CV) formula. One can 
also consider dS or AdS black holes in the bulk, such that the brane equations of motion can 
be casted in the form of Friedman equations. By introducing brane matter the cosmological 
equations can be written in the CV form which is related to AdS black hole entropy (these and 
various other aspects have been discussed in [56] quiet extensively). Now we will generalize these 
results to fijVj gravity model. 

4 Generalization to bulk f(R) gravity 

In recent years modifications of Einstein-Hilbert action by higher curvature terms is a subject 
of great interest. A very promising candidate among such modifications is the f{TZ) gravity 
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theory. The modifications introduced by the f(TZ) term in the Lagrangian can address variety 
of problems. This theory also has the potential to survive all known tests of general relativity. 
It is therefore worthwhile to explore the nature of the effective gravitational field equation on a 
brane where the bulk is endowed with f{TZ) gravity. 

4.1 Effective Field Equation on the (n-1)-dimensional Brane 

We consider a n-dimension bulk spacetime endowed with f{TZ) gravity. Following the standard 
procedure i.e., starting from the Gauss-Codazzi equation, subsequently using the Z 2 symmetry, 
we can obtain effective gravitational field equation on the brane. This has already been derived 
earlier in the context of a 4-dimensional brane embedded in a 5-dimensional bulk (see [47] and 
the references therein). The effective gravitational field equation on the (n — l)-dimensional 
brane turns out to be: 


(85) 

where the stress-energy tensor appearing on the right hand side can be decomposed into 

several small pieces. These include effective (n — l)-dimensional cosmological constant induced 
from the n-dimensional cosmological constant, brane energy momentum tensor and its higher 
order contribution. Also due to the presence of f(TZ) gravity in the bulk there is an extra term 
and all these terms along with the non-local bulk Weyl tensor turns out to be: 

^ = —^n-lhaji + Tap H- ■!p—TaP - 5 --2- QaP (86) 

^n-1 ^n-1 ^n-1 


where the additional term due to f{TZ) gravity is Q^v. This additional term has the following 
expression: 


Q 




'1 /(^_ 1 20/^ 

.4/'(7^) 4 3 /'(7^) 


2VaV^/^] , 2 VgVbf (7^) , 

3 /'(7^) J '"'" 3 /'(7^) 


(87) 


Note that in the Einstein-Hilbert limit i.e. f{TZ) TZ the above term identically vanishes and 
the effective equation reduces to Eqs. (14) and (15). In the expression of the energy momentum 
tensor, A„_i is the (n — l)-dimensional brane cosmological constant (for explicit expression see 
Eq. (16)), T^i, is the brane energy momentum tensor, iTap contain higher order terms like rj^rp^ 
etc. (see Eq. (13) for detailed expression). Finally represents the non-local bulk effect. 

In order to simplify the expression for Qfj,v, we can impose some assumptions. The simplifi¬ 
cation can be significant with the choice: dfj,TZ = 0. Then using Taylor series expansion of bulk 
curvature TZ around y = 0 hypersurface leads to: TZ = TZq + TZiy -t- 7?.2?/^/2 -|- 0{y^). As we 
have assumed that the bulk curvature depends only on the extra dimension y, all the coefficients 
appearing in the Taylor series expansion are constants. Hence we can conclude that all the 
derivatives calculated at y = 0 yield a constant contribution which does not depend on any of 
the brane coordinates. In this case we can rewrite the stress-energy tensor as: 

4 2^ 

= -^eshap + Tap + -^^Tap - 5 -( 88 ) 

^n-1 ^n-1 
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where we have an effective cosmological constant, 


■^efF — ^n—1 


F{n) 


(89) 


Thus the cosmological constant gets modified due to f{TZ) gravity in the bulk leading to an 
additional constant term F{TZ) to the induced cosmological constant A„_i. This can provide a 
possible explanation to the cosmological constant problem by fine tuning A„_i with F(TZ) term 
where the quantity F(TZ) has the following expression: 


F(7^) = 




y=o 


(90) 


Having derived the effective gravitational field equation in (n— l)-dimensional brane we will now 
write down the counterpart of this relation in (n — m)-dimensional brane as well. 


4.2 Generalization to (n-m)-dimensional Brane 

In this section we will generalize the above setup to an arbitrary lower dimensional hypersurface, 
namely, the (n — m)-dimensional brane. This can be done easily by taking a cue from the 
discussion in Einstein gravity. The effective field equation in the (n — m)-dimensional brane 
leads to the following expression: 


'i) 


(91) 


with the following expression for the energy-momentum tensor as: 
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(92) 


Here An-m is the induced brane cosmological constant, represents the brane energy 

momentum tensor and is the higher order term. All the other extra terms originate 

from the effect of higher dimensional energy momentum tensor, f{TZ) gravity and non-local 
effects from the bulk. However just as in the previous section in this case also we can simplify 
the expression further by assuming that bulk curvature depends only on the extra dimension y. 
Also if we assume that there is no energy momentum tensor on the higher dimensional branes 
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and the bulk is AdS, then Weyl tensor vanishes identically. The effective gravitational field 
equation, then takes a particularly simple form. In this case following [47], we arrive at: 

(n-m)r^ . _ „2 f \ , ^ra-m+1 

— An-mq^B + AlB + 2 ^ AB ) 

V ^n—m / 

where A„_m has to be constructed from Aeff = A„_i — FiJZ) by standard procedure illustrated 
for Einstein gravity. Thus in this case the effective equation would be identical to that obtained 
from general relativity, except for a modified cosmological constant. Also is related to 

(n — m)-dimensional gravitational constant, is the matter energy-momentum tensor on 

the {n — m)-dimensional brane, with being its higher order extension. 

This completes our discussion on effective equation. However as an illustration we will con¬ 
sider a simple application of our work in a cosmological context in the next section. 


5 Applications to f(R) gravity 

Having discussed effective field equation in f(TZ) gravity theory on a lower dimensional brane, 
alike general relativity we apply this in a static spherically symmetric and cosmological context. 
For that purpose we start with 6-dimensional bulk spacetime and try to obtain a solutions from 
the effective field equations on a 4-dimensional brane. 


5.1 Static Spherically Symmetric f(R) brane 

We start by considering the case for static, spherically symmetric spacetime with no matter 
source present. Thus brane energy momentum tensor vanishes however the electric part of the 
Weyl tensor survives. Hence the effective gravitational field equation on the four dimensional 
brane starting from the six dimensional bulk can be obtained from Eqs. (91) and (92) as. 


|‘>g^b = -4 


- 2F{TZ) -I- QAB — £aB — 


FajHCj^CB T QAB 


- -E 


(94) 


where F{TZ) is given in Eq. (90). In arriving at the above equation the cosmological constant 
has been neglected. The electric part of the Weyl tensor in such a static, spherically symmet¬ 
ric spacetime can be written explicitly using two unknown functions depending on the radial 
coordinate as, 


£ab = -1 


^4 L 


U{r) \ uaub + -^qAB 


P{r) 


tatb - -^qAB 


(95) 


with qAB being the induced metric on the t = constant surface, such that u^^qAB = 0. The five 
dimensional electric part of Weyl tensor would have similar decomposition depending on two 
separate radial functions U{r) and P{r). On projecting this five dimensional tensor Eap on the 
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four dimensional brane, the effective field equations reduce to 


Gab = --^K 

3k| 




QAB 


U{r 


UAUB + -^QAB 


U{r) 


1 


UAUB + -T^qAB + P{r) 


tatb - 2^ab 


P{r) [ tatb - -^qAB 
(96) 


Now Kn and are related through Eq. (18), such that, = (16/3 )k5(1/(t), k| = ( 6 /S)k 4 
and K 4 = SttG. Using these relations the above equation can be written in terms of the brane 
tensions, cr for the five dimensional brane, E for the four dimensional brane and the Gravitational 
constant G as, 


Gab = - (^-^F{TZ) + 47 rGE ^ qab + G{r) (^uaub + + P{r) (^tatb - 
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9-\/6 (tVSttG 


U{r) ( UAUB + ^qAB ) + P{r) { tatb - ^qAB 


Now assuming a static, spherically symmetric metric ansatz of the form, 

ds^ = + r^dn^ 

we obtain various components of the gravitational field equations to be 
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[U{r)-P{r)] 


( 101 ) 


Hence the effect of multiple extra dimensions is manifested in the lower dimensional effective 
gravitational field equations through the appearance of extra terms U{r) and P{r). Also due 
to f{TZ) gravity in the bulk there is one more additional term in the effective field equations, 
namely E(7?.). In general, these equations has to be supplemented by the conservation equation, 
which will provide a differential relation connecting U(r), P{r), U(r) and P{r) which amounts 
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to, 
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3 \ dr 
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+ t(3C^ + ^) + V = ° 


( 102 ) 

(103) 


Alike the situation in Einstein-Hilbert action all the objects on the right hand side of Eq. (99) 
depend only on the radial coordinate which can be integrated leading to the following solution 
for e~^, 


.-A 


= 1 - 


2GM Q{r) 
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47rGE 
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(104) 

(105) 


In order to obtain the solution for e’^ we need to know the solutions for the dark radiation and 
dark pressure terms. For that we will use the equations of states, 2U + P = 0 and f7 + P = 0, 
such that. 


P(r) 

Q{r) 
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Then the full static spherically symmetric solution turns out to be. 


e" = e-^ = 1 - ^ ^ 

r 2r^ 


1 \ r‘^ 

--P(P)+4^GEj - 


(106) 
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This is an exact solution to the effective field equations where no assumptions or approximations 
have been made. The effect of extra spatial dimension appears through the charge term, i.e., 
term and the term and the effect of modified gravity appears through the term. Note that 
if we assume P(P) = SttGE, then our solution will coincide with the Einstein-Hilbert action. 
Hence the addition of higher curvature terms to the bulk action can provide a potential solution 
to the cosmological constant problem. 

The static spherically symmetric solution being dS or AdS has three important aspects to 
explore, namely — (a) stability of the solutions, (b) anti-evaporation and (c) black hole thermo¬ 
dynamics. We will now comment on these features. The stability of solution can be ascertained 
by studying perturbations around it. These are again of three types, scalar, vector and tensor 
modes. These modes satisfy a set of decoupled wave equations with the following structure. 

Here □ stands for d’Alembertian operator, /(r) = e'" = e“^, $ represents the perturbation modes 
and V is an appropriate potential connected to the perturbation 4>. For stability, we require the 
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potentials to be positive. It turns out that for all choices of F{TV) the scalar and tensor modes 
have positive potentials. While for the vector mode we have stability for F{TV} > SttGE, while 
the vector modes are unstable when F{TV} < SttGE. 

The anti-evaporation effect of this solution comes into picture when we set Pq = 0 = Pq. Then 
the spherically symmetric solution represents a Schwarzschild (A)de-Sitter spacetime, which 
under proper limit can be mapped to the Nariai spacetime. A black hole in the Nariai spacetime 
has the peculiar property of increasing surface area due to quantum corrections as shown by 
Bousso and Hawking [51,52]. This phenomenon of anti-evaporation subsequently was generalized 
for Nariai black holes in the context of f{TZ) gravity [53]. Here f{TZ) gravity plays the role of 
the anomaly induced effective action leading to anti-evaporation of black holes. In our case as 
well with Pq = 0 = Pq, we have Nariai black hole as a limit. Thus from previous discussions it 
is clear that our solutions will also exhibit anti-evaporation effect. However for Pq and Pq being 
not equal to zero, our solution cannot be reduced to the Nariai form and thus in general the 
solution presented here will not exhibit anti-evaporation. 

Finally, the case of black hole thermodynamics is at our hand. The horizon associated with 
the above spacetime structure has the location r = r^, which one can determine the following 
equation. 
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The mass of the black hole (or, equivalently the internal energy) can be obtained readily from 
the above equation as. 
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The entropy can be obtained from holography, i.e., S = kBA/Ah. In this particular case of 
spherically symmetric spacetime choosing h and Boltzmann constant appropriately we obtain, 
S = Thus internal energy of the system can be obtained in terms of the entropy from 
Eq. (Ill) as, 
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Then the black hole temperature can be obtained in terms of entropy as, 
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which immediately leads to the following expression for specific heat as. 


Gv = T 



^ + I# + + {\F(Jl) - 47rGE) ^ 

(ip(P) - 47rGE) ^ ^ ^ ^ 


( 112 ) 


(113) 


(114) 


It is clear from the above expression that the specific heat shows discontinuity in its behaviour, 
thereby indicating a second order phase transition. The surface of discontinuity in the specific 
heat is given by, 


ip(P) - 47rGE 


1 ^ 9aPo 
'S^ 2S 


5Po/3 

53/2 


(115) 
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In order to understand physics behind this phase transition, it is always illuminating to discuss 
some limiting cases. For example, if we assume pure Einstein’s gravity, the specific heat can 
be given by Eq. (64), which is everywhere negative. This can also be divergent, provided the 
entropy becomes negative which cannot be satisfied in general. Thus bulk terms with positive 
dark pressure cannot lead to second order phase transition. In the other limit, where Pq = 0 = ^b) 
we get the divergence of specihc heat to correspond to the condition: S~^ = {1/2)F{TZ) — 4:TtGT,. 
Thus for F{TZ) > SttGT, we have second order phase transition. Our calculations therefore 
conhrm that Schwarzschild Anti-de Sitter solution shows second order phase transition. 

After deriving the static spherically symmetric solution along with its various properties, it 
is now time to discuss some novel features of this solution: 

• The following restrictions provided by 2U{r) + P{r) = 0 and P{r) +3U{r) = 0, shows that 
the metric elements satisfy the condition e~^ = e'^. The line element so derived depends 
explicitly on the form of U{r), U(r) and the gravity model /(P). Thus the geometry of 
the spacetime is uniquely determined only when knowledge about multiple dimensions are 
obtained through U(r) and about the alternative gravity model. 

• The solution is also stable under both scalar and tensor perturbations, since the potentials 
associated with respective wave equations are positive for scalar and tensor modes. Thus 
the black hole solution, will not change its character under scalar and tensor perturbations. 

• Finally, the specific heat associated with the black hole solution shows second order phase 
transition. This shows that the black hole transforms from thermodynamically unstable 
phase to thermodynamically stable phase. Moreover for vanishing of dark radiation we 
will obtain its mapping to Nariai black hole and therefore shows anti-evaporation. 

Having discussed the effect on black hole solution, we will now discuss the cosmological solutions 
and its implication in detail. 

5.2 Cosmological Scenario 

As in the case of general relativity, in cosmological scenario as well we can apply the effective 
field equations and obtain relevant solutions. Then we can check, where they are compatible 
with standard cosmological evolution. In the previous section we discussed vacuum solutions 
with Weyl tensor providing non-zero effect. In this we will consider exactly the opposite one, 
we will take Eab = 0 i.e. contribution from electric part of Weyl tensor is vanishing but there is 
matter in the spacetime. Then effective held equation in 4-dimension takes the following form: 

Gab = 47rGE^ qab + + f^t'^AB 

+ + QAB |ta/3s“s^ - 

+ ^4 + qab (116) 

In this expression along with brane energy momentum tensor tab and its higher order term 
two additional contribution from 5-dimensional matter energy-momentum tensor is present. 
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We assume that the brane is filled with perfect fluid with the following form for the energy- 
momentum tensor on the 4-dimension as: 

= diag(-p,p,p,p) (117) 

where p represents energy density of the matter fields and p yields its pressure. The energy- 
momentum tensor in 5-dimension is taken to be pressure free dust, such that its energy momen¬ 
tum tensor has the simple form: 


=diag(-p 5 , 0 , 0 , 0 , 0 ) (118) 

where the subscript ‘5’ indicates that it is matter from 5-dimensional brane. Thus we arrive at 
the following expression for the higher order components: 

n« = n = -^pI (119) 

Having obtained all the components in the above effective field equation, we can obtain the 
effective field equation itself. Both the time-time and spatial components of effective equation 
lead to: 

3H^ = — ^—-F{TZ) + AttGY^ + k^A4 + k\p + 

+ 2 ^ = - ^-iF(7^) -b 47rGS^ - kIp - {p^ + 2pp) - ^nlp^ - ^4 ,pI ( 120 ) 


For the bulk matter we have the following relation, namely, p 5 = {p^p /where i rep¬ 
resents the finite length of the 5-dimensional brane coordinate. After obtaining the Friedmann 
equations on the brane let us write them down in terms of relative abundances Hi = pipjPc- The 
energy density appearing on the right hand side of Friedmann equations can be expanded in terms 
of baryonic and relativistic matter fields. For higher dimensional matter being non-relativistic 
the Hubble parameter turns out to be. 


= Hi 


Ha + (Hb + H5)(^)' + H«(^)' 


( 121 ) 


where, H 5 = (K|/167rGP)(p5,o/pc) and Ha is from the combination [—{1/2)F{JZ) -|-47rGS]. Thus 
we have reproduced the standard ACDM model. At late times Ha would dominate leading 
to accelerated expansion of the universe, however this acceleration is originating since the bulk 
spacetime is being described by f{TZ) gravity theory and we are living on a brane with nonzero 
brane tension (the S term). Starting from the 1/a^ term, leading to radiation dominated uni¬ 
verse, we have matter (i.e., non-relativistic matter) dominated epoch such that ^ a~^ and 
finally the accelerating sector. The intriguing fact about the matter sector is that along with 
the baryonic matter, due to existence of extra spatial dimensions we have an additional contri¬ 
bution. This additional contribution might act as a dark matter candidate. In standard ACDM 
cosmology, the dark matter has to be included by hand while in our model it appears naturally. 
To get an estimate of five dimensional parameters we can use the epoch of matter, radiation 
equality leading to, 

1 + Zeg = 3.9 X 10^ ( 0.022 + ) (122) 

\ IoTtG-C Pc ) 
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where h = 0.72, ^ stands for the length of the extra dimension. From the redshift Zeq one can 
determine the value of the higher dimensional gravitational constant k\. On the other hand, 
unlike Einstein-Hilbert bulk in this case as well by tuning the f{TV) model itself it is possible to 
have inflation. Thus the effective field equations on the brane are capable to explain the various 
phases of the universe. 

Let us now explore the effect of additional correction terms on the Hubble parameter. To 
confine ourselves to one representative example, if we assume the matter in the 4-dimensional 
brane to be pressure-less then we have. 




1 

3 


-^F{n)+4TrGE 


8ttG po 1 4 P 0 I 1 2P5,o 1 1 iPlfi 1 

3 a3 36'^^a6 £ a3 64^® £2 ^6 


(123) 


where po represents the value of the energy density in the present epoch and p 5 _o represents the 
same for 5-dimensional matter. The above expression for Hubble parameter can be simplified 
and arranged properly leading to: 


tt2 _ ^ ( 1 W(T?\ _L A ^ 2 7*5,0 42,1 4 p|_o \ 1 

^ - “q ( + ) -7 |-^Po + r 73 + 1 + 77^6— ) 76 
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(124) 


Following the techniques developed in [49] we can introduce three constants, namely: 


8ttG 1 2 P 5 ,o, 


142,1 4^5,0 


Cf = —A -:^F(7^) +47rGS ) ; G, = —po + Ggq = -44 + 


6 “ £ 

such that the following solution to the scale factor is obtained, 

24 C] 4 Gfa^ + Gpa^ + Cs^ + 2Cfa^ + Gp 


8Q 


exp 


3\/CV (t- to) 


24CJ4Gf + Gp + Gsq + 2Cf + Gp 


Gf>0 


(125) 
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3y/—Cf {to — t)= sin ^ 


2G/C 


G„ 


Gl - 4G/Gs, 


— Sin 


2C/ + 

Gj - ACfGso 


Gf < 0, Gl > 4G/Csq 


(127) 


Though the solutions are a bit complicated, we can impose certain conditions on them such that 
these solutions simplify substantially. For that purpose we consider the situation, where G/ = 0, 
i.e., F{'R) = SttGE. Since E is positive definite the above result suggests that FijVj is positive. 
Applying this for matter dominated epoch, the Hubble parameter from Eq. (124) turns out to 
be. 




sttg r 4 P5 
3 ^ 167rG p 


47rGp / 

F{TZ) { ^ 16 4 p4 


(128) 


Assume p to be energy density of non-relativistic matter we arrive at, 
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leading to the following solution for the scale factor: 


a 


3 




U-kGpI ( 9 pi,o \ 


t 


(130) 


It is clear from the above expression that the universe undergoes a transition in the expansion 
rate. The time scale when this happened depends on whether the brane matter energy dominates 
over the bulk matter or not. Below we present the result for timescale of transition for both the 
situations: 


^ ^ \/SttGE 

» p^) ( 132 ) 

«5 V «5P5.0 J 

Thus at early universe we have a high energy regime, where a ^ while at late time low 
energy regime the scale factor variation with time modifies to a ^ which is the standard 
evolution of the matter field. Note that in Einstein gravity as well e had similar transition in the 
scale factor. But in this particular case, importantly the transition time depends crucially on the 
f{TZ) model we are considering. Similar phenomenon will be observed for radiation dominated 
epoch as well. 

The key aspects of the above results are summarized as follows: 

• First and foremost, the issue of dark energy in this particular case is resolved by the 
addition of higher curvature terms to the Einstein-Hilbert action and through the brane 
tension. The same f{TZ) model can also explain the inflationary scenario without invoking 
any scalar field in the picture. 

• The cosmology as presented by modified Friedmann equations automatically lead to an 
additional matter component originating from higher dimensions, which behave as a dark 
matter candidate. 

• Finally, note that due to presence of additional correction terms in the effective equations 
on the brane few new cosmological phases enter the picture. For example, the matter 
dominated epoch gets divided into two parts and this depends explicitly on the f{TV) 
model under consideration. This modifications will have indirect effects on the structure 
formation, which can be used to constrain not only higher dimensional parameters but also 
the form of alternative gravity theories. 

Thus with matter present in the 5-dimensional brane we can have standard cosmology with 
proper scaling of the scale factor with time. Moreover from the effective field equation it is 
clear that the 5-dimensional matter adds to the 4-dimensional one and enhances the total non- 
relativistic matter content of the universe (4-dimensional brane we live in). Structure of angular 
power spectrum, perturbation in this multiple extra dimension scenario can lead to interesting 
and important results by providing constraints on various parameters in this models. 

It is possible to ascribe the cosmological solutions so obtained holographic interpretations. 
Following [54] it suffices to show that the usual four dimensional cosmological equations can 
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be re-written as an two dimensional conformal field theory (CFT). For this we can define the 
Hawking, Bekenstein and Bekenstein-Hawking entropies as, 


{n-l)HV ^ {n-l)V 
“ 4G ’ ~ 4Ga ’ 


Sb = —£ 

n 


(133) 


where E = pV is the total energy within a volume V and pm oc but p = pm + F"(7^)/167rG. 
This is how the effect of f(TZ) gravity enters the picture, by altering the definition for energy. 
Furthermore the total entropy is by the Hubble entropy Sh, which is the entropy of the black 
hole with the radius of the Hubble size. Using these definitions the Friedman equations can be 
rewritten as. 


Sh = - VEbh (2U — KEbh)', KEbh = n{E + pV — ThSh) (134) 

n 

with energy of the black hole being Ebh = n{n — VjV/STrGa? and Hawking temperature being 
Th = —H/2ttH. This is exactly the form one would have from the CV entropy relation. It 
explicitly shows the connection between four dimensional classical gravitational physics and 2- 
dimensional quantum 2-dimensional CFT [55]. One can also consider dS or AdS black holes 
in the bulk, such that the brane equations of motion can be casted in the form of Friedman 
equations. By introducing brane matter the cosmological equations can be written in the CV 
form which is related to AdS black hole entropy [56] . 

6 Discussion 

Motivated by the success of extra dimensional models to explain various physical problems, e.g., 
hierarchy problem, cosmological constant problem, in this work we have constructed and worked 
within the premises of a very general higher dimensional model. Since gravity can propagate 
in the full bulk spacetime, the field equations for gravity on a lower dimensional brane is taken 
to be induced from the higher dimensional one. Previous works in this direction involve two 
assumptions, namely, the bulk is five-dimensional while the brane is four-dimensional and gravity 
in the bulk is originating from Einstein-Hilbert action. However in the spirit of string theory 
it is pertinent to ask what happens to this induced gravity program as the bulk spacetime is 
n-dimensional while the brane is (n — m)-dimensional, for arbitrary m. Also from the recent 
success of f{TZ) gravity it is interesting to discuss the situation when the bulk spacetime has 
f{TZ) as the action not TZ. Below we summarize all the results obtained in this work along these 
lines: 

• Starting from a n-dimensional bulk spacetime endowed with Einstein gravity, we first derive 
the effective field equations on (n — 1) and (n — 2) spacetime dimensions. Taking a cue 
from this analysis we could read off the effective field equations in (n — m)-dimensions for 
arbitrary m. 

• As an application of this result, we have illustrated two situations. One, in which there is 
no matter but bulk effects populate the right hand side of effective field equations. The 
static spherically symmetric solution so obtained has the standard 1 — (2GM/r) term, 
with additional terms originating from the bulk Weyl tensor. These additional terms 
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have contributions from both the five-dimensional and six-dimensional Weyl tensors (in 
particular its electric part). Secondly, cosmological situation with matter but no bulk 
geometric effect. From the solutions so obtained it becomes evident that matter in higher 
dimensions can add to the matter content of the visible brane and has the potential to act 
as dark matter source. More importantly, the effective Friedman equations can be casted 
in a holographic language, such that it can be related to two dimensional CFT. 

• Modifying Einstein-Hilbert action in the bulk by addition of f{TZ) term is the next im¬ 
mediate generalization. In this context as well, we have first derived the effective field 
equations in (n — l)-dimensional brane starting from n-dimensional bulk. Then we have 
obtained the effective field equations on a (n — m)-dimensional brane, starting from the 
original n-dimensional bulk. 

• Alike the situation in general relativity, in this context as well we have discussed two cases. 
For the spherically symmetric, static vacuum spacetime, f(TZ) gravity modifies the general 
relativity solution significantly. The solution also exhibits two remarkable effects. Firstly, 
it shows divergent specific heat, which points out the effect of phase transition. Secondly, 
it can be written in Nariai form, which shows the anti-evaporation effect. 

• In the cosmological context as well f{TV) gravity modifies the cosmological constant and 
the transition time scale depends explicitly on the f{TV) gravity model under considera¬ 
tion. Also the f{TZ) gravity model exhibits the CV entropy relation, such that Friedman 
equations can be rewritten in terms of two dimensional CFT. 

Thus in a nutshell, starting from Einstein’s theory in the bulk, which is taken to be n- 
dimensional we have constructed various induced objects on any arbitrary (n — m)-dimensional 
brane providing a complete generalization of existing set up. We generalize this analysis further 
by incorporating f{TZ) gravity in the bulk instead of Einstein’s gravity. In f{TZ) gravity as well 
we follow the same procedure by relating geometrical objects like extrinsic curvature, Riemann 
curvature in the n-dimensional bulk with the (n — m)-dimensional brane. In this case as well 
from the structure of the effective equation we can immediately generalize the set up on an 
arbitrary lower dimensional hypersurface as well. 

Having derived all the necessary theoretical ingredients, we consider possible applications 
of our results in both the gravity theories. As an illustration, we took a 6-dimensional bulk 
spacetime in which a 4-dimensional brane is embedded. In static spherically symmetric vac¬ 
uum spacetime we have additional terms from bulk Weyl tensor as the source of gravity. For 
f{TZ) gravity we have an equivalent of cosmological constant making the solutions asymptoti¬ 
cally (Anti) de-Sitter. In cosmological context the 5-dimensional sub-manifold contains matter 
showing that the standard cosmology to be an inherent and natural consequence of our model. 
Finally we comment on the possibility of the bulk matter as a possible candidate of dark matter 
in the lower dimensional brane. 
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A Appendix: Detailed Calculations 

In this section, we provide some detailed expressions, which have been used in the main text for 
arriving at various results. 

A.l Identities for the Derivation of (n-1)-dimensional Effective Field 
Equation 

The starting point of obtaining the effective field equation is the Gauss-Codazzi equation, which 
relates the (n — l)-dimensional curvature tensor to n-dimensional curvature tensor. For our 
purpose these two tensors are related through the following relation: 



(135) 


where we have used the definition: 


Kap = e’^epVanb 


(136) 


with e = riin®. Contraction of Eq. (135) with the (n—1) dimensional metric ha^ = e“e^ {gab — ^naUb) 
leads to the connection between (n — 1) dimensional Ricci tensor to n dimensional curvature ten¬ 
sor as, 


= ‘'^'‘Rabcdey^’pelethP'^ - e {K^^aKp^ - h'^'^ 

= ^^^Rabcdeyi (g“ - - e (iC„,7F; - 

= - e Rabcde^n^- e - Ko^^K) 


(137) 


The Ricci scalar can now be obtained as, 


= Racelelh^^^ - e Rabcdc'^.n^eln^h^^^ - e - K^^K) 

= {g^^ - enV) - e ^^^RabcdU^n^ {g^^ - enV) - e - K^) 

=(") ^_2e ^^^Racn‘^n^ - e - K^) 


(138) 


The above relation can be further simplified by using an identity for Rabn°'n^, which can be 
written as: 


^^^Rabn‘^ri^ = Va (n^V,,n“) - Van^Vbn'^ - Vb + {VbU^f 

= -V, {Kn^ - a*) - 


(139) 


Then we obtain ^'^^R in terms of ^^R as: 


(”)R R - e {Kf.aK^'^ - K^) - 2eV, {Kn^ - a*) 


(140) 
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The relation between (n — 1) dimensional Einstein tensor and n dimensional Riemann tensor is 
therefore given by 

(n-l)^ _ (n-l)D — -h ("- 1 ) 7 ? 

= - e - e - K^^K) 


2 


Mr- 2g - K^) 

= ( - \gae e“e= + £ RacU^n^K^ - £ Rabcde^n^ 


{K^ - 


(141) 


The Riemann tensor can be decomposed into Ricci tensor, Ricci scalar and Weyl curvature tensor 
in n dimension as, 


^ ^ Rated — ^ ^ Cabed T ^ {9aeRbd QadRbe QbeRad “t“ dbdRae') 
n — Z 

{gaedbd 9ad9be^ R 


(n — l)(n — 2 ) 


(142) 


which suggests, 


Rabedcln^ln'^ = ^'^^Cabedeln^e^n^ + ^ [RaeU’^n^K^ + eRaeclel) 


" h(yii,R 


(n — l)(n _ 2 )' 

Then substitution of the above equation in Eq. (141) leads to 

("-i)G„; 3= WGa 6 e“e^ + £(^) ^ 

Y Tlj Ld J Tl Zd 


(143) 




1 


(n — l)(n — 2 ) 


+ £ 
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(144) 


where we have defined, Cabedean^e^n^. This yields the n dimensional Einstein 

equation, with the following expressions: 

^^^Rab - hab R = -= kI (Tab - -^OabT) 

2 n — 2 \ n — 2 J 

(145) 
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We therefore have the following expression 
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(146) 


This is the result used in Sec. 2.1. 


A.2 Identities for the Derivation of (n-2)-dimensional Effective Field 
Equation 

We start with the connection between (n — 2) dimensional curvature tensor and (n — 1) dimen¬ 
sional curvature tensor as, 

Rabcd = RaPfiveAe-B^^D ~ {T^adK-bc — J^acK-bd) (147) 

where we have the extrinsic curvature on the {n — 2 ) dimensional surface as: K-ab = aSp, 

with Sa being the normal to the surface. Then the induced {n—2) dimensional metric turns out 
to be, qAB = e^e^(/ia /3 — SaSp). Contracting the above equation with leads to, 

- (/C^^/Cg - JCacK:) 

= {h^'' - s^s4 - {ICadJCS - ICacK) 

= - {KadK.^ - JCacIC) (148) 

On further contraction we arrive at the following result, 

0 - 2 )i? = qAB ("- 2 ) 7 ^^^ 

= - 2 - {JCabK.^^ - K .4 (149) 

Then we can use Eq. (139) to obtain an identical relation for This leads to the 

relation: 

(n-i)^ =0-2) ji _ - /C2) - 2L»„ (/Cs“ - s^Dps°‘) (150) 

which on using Eq. (140) yields: 

(”)i7 R - {KabK.^^ - T?) - ‘2Do, (/Cs“ - s^Dps°‘) 

+ - K^) + 2V, {Kn^ - a*) (151) 
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In order to obtain the effective equation on this (n — 2) dimensional surface we need various 
contractions of the Riemann tensor. For that purpose we start with the following identity, 




(n - 2 ){n - 3) “ Kvhp^) 

£aB + [qAB 


1 

(n — 2)(n 




(152) 


where we have defined the tensor Sab =^” C'a/s^ii/ej^s^e^s'^. Then the effective equation on 
the (n — 2) dimensional hypersurface turns out to be, 


(n-2)^ 


AB 


+ ^’^-^'>Ro.0s‘^s^qAB - Sab 

(,AB 3 ) '?■*'’ ‘ 

ICabIC - K.acK.% - \qAB (/C^ - ICabJC^^) 


1 

(n — 2)(n 
ICabIC- 


— ^qAB ^^R — Sab 

ICacK-% — -j^qAB — ICcdKP^) 


(153) 


Now instead of we can use the following relation: ^'>Rap = ^'^Gajs + ^hap ^'^R, 

to replace Ricci tensors in terms of Ricci scalar and Einstein tensors. This leads to, 


^^-^^Gab = 


(n-4) 


(n-3) _ 
~ Sab + 




qAB 


ICabK. — K-acJ^b ~ 2^^^ ~ ICcdKP^) 

Now introducing the Einstein equation in (n — 1) dimension we arrive at 

2 




(n-3) 


^-iT 


(154) 


(155) 


Using this expression in Eq. (154) leads to the effective field equation in (n —2)-dimension, which 
is presented in Sec. 2.2. 
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